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Abstract 

In this paper, the effect of a positive cosmological constant on 
spherically symmetric collapse with perfect fluid has been investigated. 
The matching conditions between static exterior and non-static inte- 
rior spacetimes are given in the presence of a cosmological constant. 
We also study the apparent horizons and their physical significance. It 
is concluded that the cosmological constant slows down the collapse of 
matter and hence limit the size of the black hole. This analysis gives 
the generalization of the dust case to the perfect fluid. We recover the 
results of the dust case for p = 0. 

Keywords : Gravitational Collapse, Perfect Fluid, Cosmological Con- 
stant 

1 Introduction 

The cosmological constant is an energy associated with the vacuum, i.e., 
with empty space. The inclusion of the non-zero cosmological constant into 
the Einstein field equations has been discussed several times in the past 
for theoretical and observational reasons [1]. First, it has been introduced 
by Einstein to save the universe from expanding and rejected by him after 
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expansion has been discovered by Hubble. The results of type la supernova 
[2,3] show that the universe is accelerating rather than decelerating. These 
results suggest that our universe can have a non-zero cosmological constant. 
Another analysis [4] of the peculiar motion of low-red shift galaxies give 
further evidence for the possibility of finite cosmological constant. These 
results have increased interest to study the properties of the universe with a 
non-zero cosmological constant. The physical applications of a cosmological 
constant are huge, restricting not only the growth of the universe but also 
the structure formation and age problems. The cosmological constant affects 
the properties of spacetime and matter. Since the metric of the spacetime 
and the stress-energy tensor of matter are related through the Einstein field 
equations, the effects of a cosmological constant can be analyzed by specifying 
the metric and the stress-energy tensor. We study gravitational collapse to 
see such effects. 

Gravitational collapse is one of the important issues in General Relativity. 
This theory predicts solutions with singularities and such solutions can be 
produced by the gravitational collapse of non-singular, asymptotically flat 
initial data [5-7]. Spacetime singularities can be classified into two kinds 
whether they can be observed or not. A spacetime singularity is said to be 
naked when it is observable to local or distant observer. If such singularity 
can reach the neighboring or asymptotic regions of spacetime, the singularity 
is called locally or globally naked singularity. A spacetime singularity which 
can not be observed is called a black hole. Is such a singularity formed in 
our universe? Penrose [8] proposed so-called the cosmic censorship conjec- 
ture to resolve this problem. According to this conjecture, the singularities 
that appear in the gravitational collapse are always covered by an event hori- 
zon. This conjecture has provided a strong motivation for researchers in this 
field. The compact stellar objects such as white dwarf and neutron star are 
formed by a period of gravitational collapse. It is interesting to consider the 
appropriate geometry of interior and exterior regions and determine proper 
junction conditions which allow the matching of these regions. 

Most of the problems related to gravitational collapse have been discussed 
by considering spherically symmetric system. The gravitational collapse of 
dust was first shown by Oppenheimer and Snyder [9]. They studied collapse 
by considering static Schwarzschild in the exterior and Friedman like solution 
in the interior. Many people [10-14] extended the above study of collapse by 
taking an appropriate geometry of interior and exterior regions. Markovic 
and Shapiro [15] generalized the work done by Oppenheimer and Snyder [9] 
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in the presence of positive cosmological constant. Later, Lake [16] general- 
ized the results of Markovic and Shapiro [15] for both positive and negative 
cosmological constant. Cissoko et al. [17] discussed explicitly gravitational 
dust collapse with positive cosmological constant. Recently, the same work 
has been generalized by Ghosh and Deshkar [18] for higher dimensional dust 
collapse with cosmological constant. 

In this paper, we discuss the gravitational collapse with cosmological con- 
stant for perfect fluid case. It is verified that our results reduce to the dust 
case as given by Cissoko et al. [17]. The paper is outlined as follows. In 
section 2, we give the junction conditions between a static and a non-static 
spherically symmetric spacetimes. Section 3 yields the spherically symmet- 
ric perfect fluid solution of the Einstein field equations with a cosmological 
constant. In section 4, we discuss the solution with some assumptions. Sec- 
tion 5 is devoted to investigate the apparent horizons and the role of the 
cosmological constant. Finally, we summarize the results in section 6. 



2 Junction Conditions 

We consider a timelike 3D hypersurface S, which divides 4D spacetime into 
two regions interior and exterior spacetimes, denoted by V + and V~ respec- 
tively. For the interior spacetime, we consider spherically symmetric system 
given by 

dsl = dt 2 - X 2 dr 2 - Y 2 (d6 2 + sin 2 Ode/) 2 ), (1) 

where X and Y are functions of t and r only. For the exterior spacetime, we 
take the Schwarzschild-de Sitter metric, 

ds\ = FdT 2 - ^dR 2 - R 2 (d9 2 + sin 2 6d(f) 2 ), (2) 
F 

where 

FW = 1 _ (3) 

M is a constant and A is the cosmological constant. According to the junc- 
tion conditions [19,20], it is assumed that the first and second fundamental 
forms from the interior and the exterior spacetimes are the same. These 
conditions can be expressed as 
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(i) The continuity of the first fundamental form over E gives 

(ds 2 _) s = (ds 2 + h = ds%. (4) 

(ii) The continuity of the second fundamental form over E gives 

[K ab ] = Kt b -K- b = 0, (a, b = 0,2, 3), (5) 
where K ab , the extrinsic curvature, is given by 

^ = -«J(^ + ^|i|f), (-.^- = 0,1.2,3). ( 6) 

Here the Christoffel symbols Y c are calculated from the interior or exterior 
metrics (1) or (2), are the components of outward unit normals to E in 
the coordinates x u ±. The equations of hypersurface E in the coordinates x± 
are written as 

f-(r,t) = r-r s = 0, (7) 
f + (R,T) = R-Rz(T) = 0, (8) 

where is a constant. 

Using Eq.(7) in (1), the metric on E takes the form 

(dsl ) s = dt 2 - [F(r s , t)] 2 (d6 2 + sin 2 9d<j) 2 ). (9) 

Similarly, Eqs.(2) and (8) yield 

(ds 2 + h = [F(ife) - ^y(^) 2 ]^ 2 - R*\d6 2 + sin 2 9d0 2 ), (10) 
where we assume that 

so that T is a timelike coordinate. From Eqs.(4), (9) and (10), it follows that 

ife = y(r E ,i), (12) 

^-Fk)^^^" (13> 
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Now from Eqs.(7) and (8), the outward unit normals in V and V + , respec- 
tively, are given by 

n~ = (0,X(r s ,t),0,0), (14) 

n+ = T, 0,0), (15) 

where dot means differentiation with respect to t. The components of the 
extrinsic curvature K^ b are 

K~ = 0, (16) 

YY' 

K^ 2 = csc 2 M^ = (— ) Ej (17) 

K = {R f-TR-^ + ^TR%, (18) 

K+ = csc 2 ^3 + 3 = (Fi?T) E . (19) 
The continuity of the extrinsic curvature gives 

K+ = 0, (20) 

^22 = ^2 + 2- (21) 

When we use Eqs.(16)-(21) along with Eq.(3), the junction conditions turn 
out to be 

(XY' - XY% = 0, (22) 

M = (— — —Y 3 + -Y 2 - -^Y'%. (23) 
v 2 6 2 2X 2 ; v ; 

3 Solution of the Field Equations 

The Einstein field equations for perfect fluid with cosmological constant are 
given by 

R^v = 8n[(p + p)u^u u + -(p - p)g txv \ - Ag^ u , (24) 

where p is the energy density, p is the pressure and = 5° is the four- 
velocity in co-moving coordinates. These equations for the line element (1) 
take the form 

Roo = -^-2y = 4n(p + 3p)-A, (25) 



X X Y 2 Y" X' Y' 

ft. = -x- 2 xF + ^(y-TF) = 4 ^-;»; A - < 26) 

y /^\2 xy 1 .v ,y <2 y /^"\2t 
-n.22 — —^7 — (^7 J — ~x? ^~y y~ ~~ ~~x ~y ~ y 

= 47r(p-p)-A, (27) 
^33 = sm 2 9R 22 , (28) 

y' y' 

«oi = - 2 y + 2 xY = °- ( 29 ) 

Now we solve these equations. When we integrate Eq.(29) w.r.t. t, we get 

X = £, (30) 

where W = W(r) is an arbitrary function of r . Using this value of X in 
Eqs.(25)-(27), it follows that 

y y n 1 - vu 2 

2 7 + ( 7 ) 2 + ^- = A-87,p. (31) 
Integrating this equation w.r.t. t, it turns out that 

m y^ 

y2 = W ,2_ 1+2 _ +(A _ 87rJo) _ ; (32) 

where m = m(r) is an arbitrary function of r and is related to the mass of 
the collapsing system. When we use Eqs.(30) and (32) in (25), we obtain 

m! = 47r(p + p)Y 2 Y' - -8np'Y 3 . (33) 

3 

For physical reasons, we assume that density and pressure are non-negative. 
Integration of Eq.(33) w.r.t r yields 

m(r) = 4tt f (p + p)Y 2 Y'dr - — f p'Y 3 dr. (34) 
jo 3 jo 

Here we take constant of integration to be zero. The function m(r) must 
be positive, because m(r) < implies negative mass which is not physical. 
Using Eqs.(30) and (32) into the junction condition (23), we obtain 

M = m- ^Y 3 . (35) 
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We see from Eq.(3) that the exterior spacetime becomes the Schwarzschild 
spacetime for A = and M as the total energy inside the surface S due to its 
Newtonian asymptotic behaviour. The total energy M(r, t) up to a radius r 
at time t inside the hypersurface E can be evaluated by using the definition 
of the mass function [10,17]. For the metric (1), it takes the following form 

M(r, t) = ir 3 i? 23 23 = \y[1 - (^) 2 + y 2 ]. (36) 
Using Eqs.(30) and (32) in Eq.(36), it follows that 

Y 3 

M(r, t) = m(r) + (A - 8?rp) — . (37) 

4 Solution With W(r) = 1 

In this section, we consider the case A — 8ivp > and the condition 

W(r) = 1. (38) 

Using Eqs.(30), (32) and (38), we obtain the analytic solutions in closed form 

as 

Y(r,t) = ( A _^ )*sinh§aM), (39) 

j . . 6m . i rr m' IQnmp' , . , 

X(r,t) = (-. — ) 5 {^— + TT r-^ismha^t 

v ' ; v A-8ttj9 ; 11 3m (A-87rp) 2j v ' 

, — 8irp'(to — t) , A — 8np^ , . . , -i . . 
+ { 1 y =L + to J -} cosh a(r, t)] smh~ a(r, t), 



y/3(A - 8np) 



(40) 



where 



J3(A - 8itp) 

<*(r,t) = V 2 " Ur)-t\. (41) 

Here to(r) is an arbitrary function of r. In the limit p — > the above 
solution corresponds to the Tolman-Bondi solution [21] 

limFM) = tot,,-*)']*, (42) 

lim X(r,t) = — — r-^. 43 

P ^ [6m^ -# 
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5 Apparent Horizons 



When the boundary of trapped two spheres is formed, we obtain the appar- 
ent horizon. Here we find this boundary of the trapped two spheres whose 
outward normals are null. For Eq.(l), this is given as follows 

gTY^Yr = Y 2 - (^) 2 = 0. (44) 

Using Eqs.(30) and (32) in Eq.(44), we obtain 

(A - 8np)Y 3 - 3Y + 6m = 0. (45) 

The solutions of the above equation for Y give the apparent horizons. For 
A = 8np, it becomes the Schwarzschild horizon, i.e., Y = 2m. When m = 
0, p — 0, it yields the de-Sitter horizon Y = The case 3m < ^=rg^; 

leads to two horizons 

2 ib 

r > = 7^w P ox r (46) 
y * = -7dW cos f-^ sin f>> (47) 



cos?/> = -3m^A - 8np. (48) 



where 



For m = 0, it follows from Eqs.(46) and (47) that Y\ = y A _ 8np and Y 2 = 0. 
Y\ is called the cosmological horizon and Yi is referred to be black hole 
horizon which can be generalized for m ^ and A ^ 8ixp respectively [22]. 
It is mentioned here that both horizons coincide for 3m = , K 1 — , i.e., 

Y 1 = Y 2 = ■ = 1 = Y (49) 

which gives a single horizon. It is obvious that the range for the cosmological 
horizon and the black hole horizon turns out to be 



< Y 2 < 1 < Y 1 < J - 3 Q . (50) 



VA - 8vrp ~ ~ V A - 8np 

its largest proper area 4ttY 2 
mological horizon has an area between . A l to . 12 J r . For 3m > n = — 



The black hole horizon has its largest proper area AtcY 2 = ^ and the cos- 
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there are no horizons. The formation time of the apparent horizon can be 
calculated with the help of Eqs.(38), (39) and (45) and is given by 

sinlr 1 ^ - 1)3, (n=l,2). (51) 



^/3(A - 8np) ' 2m 



In the limit p — > we obtain the result corresponding to Tolman-Bondi 
[21] 

4 . . 

t a h = *o - (52) 

From Eq.(51), it can be seen that both the black hole horizon and the cos- 
mological horizon form earlier than the singularity t = t . From Eq.(51), it 
follows that 

Y 

— — = cosh 2 a n . (53) 
2m 

Eq.(50) yields that Y 1 > Y 2 and also Eq.(51) gives t\ < t 2 , i.e., cosmological 
horizon forms earlier than the formation of the black hole horizon. To see 
the time difference between the formation of the cosmological horizon and 
singularity and the formation of the black hole horizon and singularity re- 
spectively, using Eqs.(46)-(48), we need to calculate the following quantities 

= V^ + ^l)< , (54) 
dm m sinip cosip 

^jJ = l(_ Sm 3 + ^ CQS 3 )>Q (55) 
dm m sm^) costp 

We define the time difference between the formation of singularity and the 
apparent horizon, denoted by r as follows 

r n = t - t n . (56) 

It follows from Eq.(53) that 

dTn 1 (57) 



d( 2^ ) sinh a n cosh a n ^3 (A — 8irp) 
Using Eqs.(54) and (57), it turns out that 

dn = _dr^d(^)_ = 1 ( _sin| + 3cos| < Q 

dm d(^) dm my^A - 8np) sinh a x cosh a x sm V> cos ^ 

(58) 
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This shows that T\ decreases and hence the time difference between the for- 
mation of singularity and cosmological horizon decreases. From Eqs.(55) and 
(57), it follows that 



ar-y 1 , sin ^ — '- 3 cos ^ — '- 



dm m v /3(A - 8irp) sinh a 2 cosh a 2 sin ^ cos^ 



) > 0. (59) 



This implies that t 2 increases which means that the time difference between 
the formation of singularity and black hole horizon increases. 



6 Conclusion 

In this paper, we have studied the gravitational collapse of a perfect fluid 
in the presence of a cosmological constant. The effects of the cosmological 
constant on gravitational collapse have been discussed in the following two 
ways. 

Firstly, the cosmological constant plays the role of repulsive force, i.e., 
it slows down the collapsing process. The cosmological term behaves like a 
Newtonian potential given by <ft = |(1 — (/oo)- Using Eqs.(12) and (35) for 
the exterior metric, the Newtonian potential takes the following form 

0( J R) = ^ + (A-87rp)^. (60) 

The corresponding Newtonian force turns out to be 

F = -^ + (A- 87 rp)| (61) 

which vanishes for R = n 1 t . and m = „ ,} „ . Thus the force becomes re- 

VA— 8irp 3vA— 8/Tp 

pulsive/attarctive for larger /smaller mass and radius respectively than these 
values. This means that the size of the black hole can be visualized by com- 
paring the repulsive and attractive forces. The repulsive force generates from 
the cosmological constant for A > 8np. It is worth mentioning here that for 
the perfect fluid A can play the role of a repulsive force only for A > 8irp 
while in the dust case this is true for all values of A > 0. From Eq.(32), the 
rate of collapse turns out to be 

777/ Y 

Y = -yi + (A-8np)j. (62) 
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For collapsing process, the force should be attractive, i.e., the acceleration 
should be negative which implies that Y < ( A 3 ^ ) 5 . Thus Eq.(62) shows 
that the cosmological constant slows down the collapsing process if A > 8irp. 
This means that, for p > the force becomes attractive and hence the 
cosmological constant does not slow down the collapsing process. 

Secondly, there are two physical horizons instead of one due to the pres- 
ence of the term A — 8irp, i.e., the black hole horizon and the cosmological 
horizon respectively. The cosmological constant influences the time differ- 
ence between the formation of the apparent horizon and singularity. We find 
that the cosmological constant affects the process of collapse and hence it 
limits the size of the black hole. In perfect fluid case, these results are valid 
only for A > 8np while in the dust case these are valid for all A > 0. Thus 
we conclude that the pressure term creates a bound for the cosmological con- 
stant to act as a repulsive force. It is mentioned here that if we take p — 0, 
the results reduce to the dust case [17]. 
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